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Abstract. Relativistic and QED corrections to the recently discovered first 
vibrational 2pau state are presented. This state has an extremely small 
nonrelativistic binding energy i?3 = 1.085045252(1) X 10~^ a.u. Its wave functions 
has a maximum at H ~ 100 a.u. and extends up to several hundreds. It is shown 
that this state does not disappear if higher order relativistic and QED corrections, 
including the Casimir-Polder effect, are taken into account. 



PACS numbers: 31.30.Jv, 31.15.Pf, 31.15.Ar 

1. Introduction 

We have recently reported |21 El on the possible existence of a new bound state of 
the molecular ion. It corresponds to the first excitation of the 2p(7„ vibrational 
state, i.e. a state with total three-body orbital angular momentum L ^ and 
with proton-proton spin coupled to S'„ = 1. Its binding energy is extremely small, 
Eb = 1.085045252(1) x 10^^ a.u. (atomic units) with respect to the p-H dissociation 
threshold, and its wave function extends over several hundreds a.u. Although not 
being populated, it results into a huge p-H scattering length value of at ~ 750 a.u., 
which dominates the low energy p-H scattering. Apart from its exotic character, the 
very existence of such a state can dramatically modify the p + 11 + hv reaction 

rate and consequently can explain the II2 molecular abundance jjfi |^ . Our previous 
calculations were done in the framework of nonrelativistic quantum mechanics. Given 
the smallness of its binding energy, it is necessary to undertake a systematic study of 
the relativistic corrections in order to see to what extent they modify, or even question, 
the results previously found. This is the aim of the present paper. 

At a first glance it may be expected that at very large distances a change of the 
asymptotic potential due to the Casimir-Polder (or retardation) effect can lead to 
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disappearance of this weakly bound state. For the charge-dipole case the asymptotic 
potential at distances r 3> 100 a.u. has a form |S] 

where is the electric dipole polarizability of an atom (for the hydrogen atom 
ad — 9/2) and a w 1/137 is the fine structure constant. The first term is the 
classical polarization potential, which results from the instant Coulomb interaction 
and is already included in the nonrelativistic solution. The second term represents 
the two transverse photon exchange between a neutral system (hydrogen atom) and a 
distant charged particle (proton) . In order to evaluate the influence of this contribution 
on the binding energy one may use the value of the retardation potential — second 
term in Eq. — at the maximum of the bound state wave function, Rmax ~ 100 
a.u., that gives: 

Kot(r™a.) = -(9/4)(lla/2^)(me/mp)i?-L ~ 10~^' a.u.. 

Thus the change of the 2pau{v = 1) binding energy is of the order Ai?i-et = Eb x 10~^. 
The smallness of this quantity suggests that relativistic and QED corrections can be 
treated using the standard perturbation technique for bound states, and eventually at 
some stage should take in the Casimir-Polder effect. 

2. Variational calculation of the nonrelativistic solution 

The numerical calculations of the bound state wave function have been performed 
using the variational approach described in details in |Zj . The variational wave function 
for an S state has the form, 

00 

*(ri,r2) = y |[/,Re[e-"''-i-'^'''^-^'''] 

^ (2) 
lni[e-"'''^-'''''^-'^'''] } - (1 ^ 2), 

where ri and r2 are the position vectors of the electron with respect to two protons. 
Complex parameters ai. Pi and 7^ are generated in a quasi-random manner: 

{A2-Ai) + Ai 

iA',-A[) + A[ , 

[a; J designates the fractional part of x, Pa and Qa are some prime numbers, [^1,^2] 
and [^2,^2] are real variational intervals which need to be optimized. Parameters Pi 
and 7i are obtained in a similar way. 

In order to get the accurate result we use 5 sets of the basis functions of the 
type (O, for which the variational parameters have been searched independently. 
The proton-to-electron mass ratio, mp — 1836.152701 rrie, has been adopted for this 
calculations. 

In Table ^ present the convergence of the computed binding energies as a 
function of N, number of the basis functions. One can see that the nonrelativistic 
binding energy for this weakly bound state has a relative accuracy of 10~^, what is 
compliant with the requirements of the precise spectroscopy. The next question is how 
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Number of state {N) 


Binding energy Eb (in a.u.) 


2000 


1.08504520 X 10"^ 


2500 


1.085045237 x lO^^ 


3000 


1.0850452464 x 10"^ 


3500 


1.0850452494 x 10"^ 


oo 


1.085045252(1) x 10"^ 



Table 1. Convergence of the binding energy (in a.u.) for the 2pc!y.(v = 1) state 
with respect to a number of basis functions 

to improve this value by taking into account the corrections imposed by a relativistic 
theory and QED. These aspects will be discussed in the following sections. As it 
was demonstrated in the introduction, these corrections can be evaluated using the 
standard perturbation expansion over the parameter a ~ 1/137, which can be derived 
from the nonrelativistic QED effective field theory 

3. Corrections due to the Breit Pauli Hamiltonian 

The Breit-Pauli Hamiltonian provide us with the relative order corrections with 
respect to the nonrelativistic energy of a state. 

The major contribution comes from the relativistic correction for the bound 
electron. 

The other term of the Breit-Pauli Hamiltonian which has to be considered, is the 
transverse photon exchange contribution, which reads: 

.p(2) / PePi , ri(ripe)pi , \ 

(4) 

I P1P2 R(Rpi)p2 \ 
2Mp2 \ R i?3 / ■ 

The remaining recoil corrections are negligibly small compared to uncertainty in the 
relativistic correction for the bound electron. 

Beyond these terms, we have included as well the correction due to the finite size 
of the proton, 

^i^nuc = ^^^^(^(rt)+^(r2)), (5) 

where Rp = 0.862(12) fm is the proton root-mean-square radius. 

As can be seen in the Table |2 the relativistic correction to the binding energy 
is of the order ~ 10"'^, what is in agreement with the work of Howells and Kennedy 
0. These authors studied the relativistic corrections for the high vibrational states 
of the Isag series in HJ. They found that the Breit-Pauli relative contribution to the 
binding energy of the weakly bound states is of the order 10~^, while in the case of low 
vibrational states it constitutes only a ~ 10~^ part of the binding energy. The value 
of AEnr in Table 13 agrees well with our previous estimate 0, based on the simplified 
approach of [5]. 
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AEnr 


-1.085 045 252 x 10^ 




AErc 


0.003 285 2(4) x 10 


-9 


AEtr-ph 


0.000 013 371 X 10- 


9 


AEfiuc 


-0.000 000 067 X 10- 


9 


AEa2 


0.003 298 5(4) x 10 


-9 



Table 2. The Breit-Pauli Hamiltonian corrections (in a.u.) to the binding energy 
of the 2pcru{v = 1) state 



The uncertainty in the relativistic correction for the bound electron is considerably 
larger than other uncertainties in corrections. This is due to the strong cancelation 
between the correction terms for the molecular ion and the ground state of the 
hydrogen atom. The very accurate variational solution, providing the accuracy for the 
nonrelativistic energy to be ~ 10~^^ a.u., is still not enough to get the precise value 
for this relativistic contribution. 



4. Radiative and higher order relativistic corrections 



The complete spin-independent contribution of order and a^{m/M) has the form 

iHnini, 

[4 



1=1,2 
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^P{L,v) + l-\]{5{v,)) 
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Ml, 



-lna-4/?(i,i;)- 
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y 



14 
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where 



{Pe{Ho-Eo) In {{Ho-Eo)/Roo)Pe 



(6) 



(7) 



47r (5(ri) + J(r2)) 

is the Bethe logarithm, Hq is the three-body nonrelativistic Hamiltonian, Pe is the 
electron momentum operator and Q{r) is the Q-term introduced by Araki and Sucher 

E], 



Q{r) 



lim 



e(r - p) 

Aur^ 



(lnp + 7£;)J(r) 



In calculating the 2pcr„(u — 1) state of H2 ^ the Bethe logarithm was taken equal 
to the hydrogenic limit, namely, (3{^S, 1) ~ 2.9841. This is justified since the electronic 
wave function for the 2pau{v = 1) state to a good extent can be approximated by the 
antisymmetrized hydrogenic wave function: V'e('"i, '"2; -R) = (l/'\/2)(V'H(ri) — ■0H(r2)))- 
On the other hand, this accuracy is sufficient to get a relevant estimate of the order 
radiative correction. 

Our calculations include also the order corrections in a non-recoil limit. Among 
them are the one-loop self-energy and vacuum polarization corrections for the bound 
electron (next to the leading term in a expansion of the external field approximation 

[12 ES]) 
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AEnr 


-1.085 045 252(1) x lO-'-* 


AE^2 


0.003 298 5(4) x IQ-^ 


AE^s 


-0.000 470 02(1) X 10-9 


AE^. 


-0.000 003 29 X 10-3 


Eb 


1.082 219 8(4) X lO"'^ 



Table 3. Relativistic and QED corrections to the 2pau{v = 1) state of the 
hydrogen molecular ion HJ. 



and two-loop QED corrections, 

se'-''^ - — 



4358 IOtt^ 



1296 



27 



The last contribution is the relativistic corrections of order a'', 6Eit\ for the 
bound electron (the ma^ order term in the expansion of the Dirac energy for the 
two-center problem). 



where Hb is the Breit-Pauli Hamiltonian for the bound electron of the two center 
problem and Q is a projector orthogonal to the initial 2pcr„ electronic state, and 



^(4) 



a 



Pe 



16m^ 
37r 
~16m4 



Pe 



i=l,2 ^ 



Pe, — 

ri 



1 

Pe, — 

ri 



128m4 



Pe 



2 1 

Pe, 



I ' fact, the relativistic correction turns out to be negligibly small and can be 
omitted. 

The summary of the relativistic and QED contributions up to and including term 
of order a* is presented in Table|31 The uncertainty in the final value is determined by 
the uncertainty in calculating the leading relativistic correction for the bound electron. 
The other corrections have been obtained with much better accuracy. 



5. Spin effects 

In the preceding evaluations, the spin effects were ignored. For the 2pau{v = 1) state, 
the spin-spin interaction has a form: 



rrHFS 



a — MeSe 



where = (1 + ae)/me and fip — {I + ap)/Mp are the magnetic moments of an 
electron and a proton, respectively. For the lowest hyperfine state, ^tot = 1/2, the 
spin-dependent correction to the binding energy is 

2 S-TT/Ze/Xp 



HPS 



5E-i 



(S„s,) (<5(ri) + <5(r2))H+ 



where S„ = Sp-^ +Sp2 is a total spin of protons, and (S„Se) = —1. On the other hand, 
the p + H asymptotic states should be antisymmetric with respect to exchange of two 
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protons as well. For the 2j3(t„(u — 1) this can be realized only when proton spins are 
parallel {Sn = !)■ In this case the threshold energy should be 



1/2) 



Then the change of the binding energy due to the spin-spin interaction is 



HFS 

+ 



HFS \ 



--1/2 



= -8.223 X 10"^"* a.u. 



6. Conclusions 

The relativistic and QED corrections to the 2pau{v = 1) vibrational state of the 
molecular ion have been evaluated. Calculations include up to order terms. 

The main conclusion of this work is that, despite the smallness of its binding 
energy, the existence of this state is not questioned by the impact of the relativistic 
and radiative effects. By including all these corrections one gets a binding energy 
Eb = 1.0822198(4) x lO'^ a.u., or Eb = 2.944 870(1) x 10"^ eV, what represents a 
relative modification AEb/Eb = 2.61 x 10~^ of the nonrelativistic value. 

While the nonrelativistic binding energy is now known to a relative precision of 
about 10~^, the final value for the physical binding energy has an uncertainty by two 
orders of magnitude larger. Generally, for the low vibrational "gerade" states of the 
hydrogen isotope molecular ions, the convergence with respect to a is better. The 
most conceivable explanation is that, for the weakly bound states, the cancelation of 
the different correction terms in the molecular ion and the atom has the strongest 
effect and this slows down the convergence of the expansion. It manifests already in 
the leading order corrections of the Breit-Pauli Hamiltonian and was observed before 
in [HI for the IstJg series of states. 

The contribution of the spin-spin interaction is rather small: its magnitude is by 
three orders smaller than the relativistic correction for the bound electron. 

Finally, we would like to mention that the existence of this weakly bound 
molecular ion state is of fundamental importance. It manifests itself in a huge p-H 
scattering length a « 750 a.u. Pj, which determines the low energy scattering of 
proton by atomic hydrogen. The formation rate is substantially influenced by the 
p-H resonant cross section. This can help to explain the abnormal abundance of H2 
molecules in the interstelar space 0] |S] . Experimental confirmation of the considered 
state would be very appreciated, despite being difficult to realize [J. 
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